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Abstract 

What is the minimum perimeter of a convex lattice n-gon? This ques- 
tion was answered by Jarnik in 1926. We solve the same question in the 
case when perimeter is measured by a (not necessarily symmetric) norm. 
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1 Introduction 

What is the minimal perimeter L n that a convex lattice polygon with n vertices 
can have? In 1926 Jarnik [I] proved that L n = ^^-n 3 / 3 + 0(n 3 / 4 ). The aim of 
this paper is to extend this result to all, not necessarily symmetric, norms in the 
plane. As usual, such a norm is defined by a convex compact set D C R 2 with 
G int D, and the norm of x G R 2 is 

IMI = \\ x \\d = mm {^ > : x G tD}. 

Let Z 2 be the lattice of integer points in R 2 , and write V n (n > 3) for the set 
of all convex lattice n-gons in R 2 , that is, P G V n if P = conv {zi, . . . , z n } where 
Zi,...,z n G Z 2 are the vertices, in anticlockwise order, of P. The D-perimeter 
of P is defined by 

n 

Per P = Per D P = ^2\\z i+1 - Zi\\ D 

i=l 

where z n+ \ = z\ by convention. Note that for a non-symmetric D, Per^P depends 
on the orientation of P as well. Define now 

L n = L n (D) = mm{Pei D P : P eV n } (1.1) 

Since D will be kept fixed throughout, we will often write PerP and L n instead 
of Per D P and L JD). 
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In this paper we determine the asymptotic behaviour of L n (D) for all norms. 
We will also show that, after suitable scaling, the minimizing polygons have a 
limiting shape. The same results were proved by Maria Prodromou [5] in 2005 in 
the case when D is symmetric, that is, D = —D. We will see that most of the 
difficulties in the non-symmetric case do not come up in the symmetric one. 

Define T as the set of all positive continuous functions r : [0, 2tt] — ► R + with 
r(0) = r(2n). Such a function is the radial function of a starshaped set in R 2 ; 
such a set contains the origin in its interior and the half-line starting at the origin 
in direction u(t) = (cost, sin t) intersects its boundary at a single point which is 
at distance r(t) from the origin. We write S for the set of all starshaped sets in 
R 2 . Every convex compact set K C R 2 with G int K is, of course, starshaped. 
We denote by T c the set of radial functions of all such convex compact sets. 

Let tq G T c be the radial function of D. The problem of determining L n (D) 
is closely related to the following variational problem, to be denoted by VP(r ). 
We seek a radial function re J 7 that minimizes 

r r 3 {t)/r (t)dt 
Jo 

subject to / r 3 (t) cos tdt = 0, / r 3 (t) sin tdt = 0, 
Jo Jo 

and - / " r 2 (t)dt = 1. (1.2) 
2 jo 

Assume r(t) is the radial function of a convex (or starshaped) compact set 
K C R 2 . Then the first condition says that the centre of gravity, g{K), of K is 
at the origin, and the second condition says that Area-fT = 1. We will explain 
later the meaning of the function to be minimized. Using the results concerning 
L n we will prove the following. 

Theorem 1.1 There is a unique solution r G T to the variational problem. It is 
the radial function of a convex compact set in R 2 defined as the only function of 
the form 

la, 

- = ho cos t + c sin t 

r r 

with a > 0, 6, c G R, that satisfies the constraints of VP(ro). 

Notice that all the positive functions of the form ^ + 6 cost + csint are 
radial functions of a convex set. Indeed, the sign of the curvature is given, in 
the differentiable case, by the sign of (-)" + - which happens to be equal to 
a((^)" + which is always positive because D is convex. This result can easily 
be extended to the non differentiable case. 

We mention further that the solution to VP{tq) is unique in a larger class 
than T . This will be clear from the proof. 



I. Bar any, N. Enriquez - Jarmk's convex lattice n-gon 



3 



2 Results and notations 

Assume that the vertices of a minimizer P n G V n are Z\, . . . , z n in anticlockwise 
order (which is the orientation giving the minimal D-perimeter) . Then E n = 
{z 2 — z\, . . . , z n — z n -i, z\ — z n } is the edge set of P n . Define C n = conv E n . Note 
that the E n determines P n uniquely (up to translation). Even more generally, 
the following is true. 

Proposition 2.1 Suppose V C R 2 is a finite set of vectors whose sum is zero. 
Assume further that u, v G V, u — Xv with A > implies that u = v. Then there 
is a unique (up to translation) convex polygon whose edge set is equal to V . 

Proof. This is very simple. One has to order (cyclically) the vectors in V 
by increasing slope as V\, . . . , v m v\. Then the polygonal path through the points 
0, v i, v i + V2, v\ + V2 + t>3, . . . , v\ + . . . + v n = in this order is a convex polygon 
with edge set V. Uniqueness is clear. □ 

We call this construction the increasing slope construction. Here come our 
main results. We let K, denote the family of all convex compact sets in R 2 with 
non-empty interior. For K, L G /C, dist (K, L) denotes their Hausdorff distance. 

Theorem 2.2 There is a unique C G K, such that limdist ((AreaC ra )~ 1//2 C n , C) = 
0. Moreover, g(C) = and \imn~ 3 ^ 2 L n (D) exists and equals 



We will prove the uniqueness part of Theorem 11.11 by showing that the radial 
function of C is the unique solution to the variational problem VP{tq). 

Theorem 2.3 There is a convex set P C R 2 such that the following holds. Let 
P n be an arbitrary sequence of minimizers, of L n (D), translated so that min{a; : 
(x,y) G P„} is reached at the origin. Then limdist (n" 3 / 2 P n , P) = 0. 

We explain in Section [TU1 how and why P is determined uniquely by C. More- 
over, it is shown in section 11 that the round shape found for P in Jarnik's case 
is obtained if and only if the unit ball D is given by an ellipse having a focus 
point at the origin. 

To avoid some trivial complications in the proofs we assume that D is strictly 
convex. We emphasize however that the above results are valid without this extra 
condition. We make another simplifying assumption, namely, that 



This is just a convenient scaling of the unit ball which leaves the set of minimizers, 
and the corresponding E n , C n and consequently C, P unchanged. 

The strategy of proof of the key Theorem 2.2 is as follows. We put together 
the following ingredients : 




Area D = 1 



(2.1) 
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• almost all primitive vectors of C n belong to E n (Section 7), 

• the normalized convex hulls (AreaC n )~ 1 / 2 C n are sandwiched between two 
fixed Euclidean balls (Section 6), so that the Blaschke selection theorem 
applies (Section 9), 

• the radial functions of the only possible limiting points of the sequence 
(AreaC ra ) _1 / 2 C n are solutions of VP(r ) (Section 5). Moreover, the varia- 
tional problem VP(r ) has a unique solution (Section 8). 

3 Auxiliary lemmas 

We write P for the set of primitive vectors in Z 2 , i.e., z = (x,y) e Z 2 (z / 0) is 
in P if x and y are relatively prime. The following two claims are very simple. 

Claim 3.1 For all n > 3, L n < L n+1 . 

Proof. Let P n+ \ = conv {zo, Zi, . . . , z n } be a minimizer for L n+i and set 
P* = conv { Zl , . . . , z n }. Then L n < Per P* < L n+1 . □ 

Claim 3.2 E n C P. 

Proof. Assume P n is a minimizer and the edge z 2 — Z\ ^ P, say. Then 
the segment [-21,-22] contains an integer z e Z 2 distinct from -21,-22- The convex 
lattice n-gon conv {z±, z, z$, . . . , z n } has shorter D-perimeter than P n because the 
triangle conv {z±, z 2 , z 3 } contains the triangle conv {zi, z, z 3 } so the latter has 
shorter D-perimeter. □ 

The following lemma will be useful when proving that most points in C n H P 
belong to E n . 

Lemma 3.3 Assume a, b G E n and a 7^ ±b. Let T be the parallelogram with 
vertices 0, a,b,a + b. If x,y e (T n P) \ E n and x 7^ y, then x + y T. 

Proof. If x + y E T were the case, then set E* = E n U {x, y, z} \ {a, b} 
where z = a + b — x — y. The increasing slope construction works now because 
J2zeE* -2 = and gives rise to convex lattice (n + l)-gon P if there is no u G E n 
with u = \z with A > 0. If there is such a u, we replace u and z by u + z in E* , 
and the increasing slope construction gives a convex lattice n-gon P. We claim 
that P has shorter Z)-perimeter than P n . This clearly finishes the proof. 

To prove PerP < PerP n we have to show that ||x|| + ||y|| + ||-2|| < ||a|| + ||6||. 
Assume that the anticlockwise angle from a to b is smaller than n. Then x,y,z G 
pos {a, b} where pos {a, b} is the cone hull of a and b. Order the vectors a, b, x, y, z 
by anticlockwise increasing slope. The outcome is a, x, z, y, b say. Then the 
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triangle A = conv{0,a, a + b} contains the quadrilateral Q = conv{0, x, x + 
z, x + y + z} so the latter has shorter D-perimeter. Now a+b=x+y+z and 



PerQ 



\x\ 



\z\\ + \\x + y + z\\ < Per A 



|a + 6||, 



and Per P < Per P n follows. 



□ 




Figure 1. The proof of Lemma 3.3 

We write B for the Euclidean unit ball in R 2 and \x\ for the Euclidean norm 
of x G R 2 . Since D is compact convex and G int D, there are positive constants 
d\, c?2 such that d\B C D C d 2 B, or, equivalently, 

d±\x\ < \\x\\ < d,2\x\, for every x G D. 

In what follows c, ci,C2, .. denote positive constants independent of n. We will 
also use Vinogradov's convenient notation: /(n) gf(n) means that there are 
positive constants c and n such that cf(n) < g(n) for all n > n . Of course, 
the constants do not depend on n. But they depend on D, more precisely, they 
depend on the constants di,d 2 . f(n) 3> g(n) has the same meaning but with 
f(n) > cg(n). We will also use the big Oh and little oh notation. 

We need some standard estimates on the distribution of lattice points and 
primitive points in a convex body K G JC, see [3] or pQ for a proof. Let L denote 
the Euclidean perimeter of K. We assume that L > 3, say, but we think of K as 
"large". In fact, in most applications L tends to infinity. The following estimate 
is simple and well-known. 



KHZ 2 \ -Are&K 



< 2L. 



This implies, with the standard method using the Mobius function, that 
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KnP\ 



-Area K 



71- 



< 3L\orL. 



(3.1) 



(3.2) 
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Assume next that / : R 2 — > R is a 1-homogeneous function, that is, /(Ax) = 
Xf(x) for every x G R 2 and A > 0. Writing M = max{|/(z)| : z G K} the 
following estimates hold. 



E /(*) - i f(z)dz 



< 2ML, 



(3.3) 



E /(*) - - 2 I 



zeA'nP 



< 3ML logL. 



(3.4) 



The same estimates hold when K is a (non-convex but) starshaped set whose 
boundary consists of finitely many line segments. (Then, of course, the perimeter 
of if is a finite number L.) This fact will be needed in Section [51 

These estimates will be used quite often in the case when K = XKq, and 
A -> oo with K fixed. Then formulae (ETI) . (ET211 . (IO> have the following 

simpler form: 

| AT n Z 2 | = A 2 Area A (l + 0(A -1 )), 
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iTnP| = - T A 2 AreaA ((l + 0(A _1 log A)). 



7T 



E /(*) 

zei^nz 2 



A 3 / /(z)^ + 0(A 2 ), 

«0 



(3.5) 
(3.6) 

(3.7) 



£ f{z) = - 2 X i j f(z)dz + 0(X 2 log X). (3.8) 

The constant in the big Oh notation depends only on K . Here K is either a 
convex set or a starshaped set with boundary consisting of finitely many line 
segments. 



4 Bounding L n 

In this section we give upper and lower bounds on L n . 

Claim 4.1 L n > n 3/2 . 

Proof. Here we use the following density principle. The sum of the lengths 
of n distinct primitive vectors is at least as large as the sum of the lengths of the 
n shortest (distinct) primitive vectors. We will see the same principle in action 
a few more times. 

Let vi , . . . , v n be the n shortest (in .D-norm) vectors in P (ties broken arbitrar- 
ily). Set A = max{||f j|| : i = 1, . . . , n}. Then (int AD) fl P C {v i, . . . , v n } C XD. 
The boundary of XD contains at most Per^A-D < 2ird2X lattice points. So 
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\\D n P| — 2nd 2 \ <n< \XD n P|. Using (JSU) with \D (recalling Are&D = 1) 
gives 

IAD n PI = — A 2 (l + 0(A _1 logA)). 

This shows that n = ^A 2 (l + 0(A _1 logA) implying that A = + o{l))n 1 / 2 . 
Using this in (13. 8p with XD gives 

n 

L n > £NI> E IN 

1 zeint (AD)nP 

4t - 0(A _1 log A)) A 3 / \\z\\dz^>n 3/2 . 

TX l J JD 



> 



□ 

Claim 4.2 L n < n 3 / 2 . 

Proof. Again, let ui, . . . , v n be the n shortest (in D-norm) vectors in P and 
set v o = — J2\ v i- By the increasing slope construction the vectors v , vi, . . . , v n 
form the edge set of a unique (up to translation) convex lattice n-gon or n + 1- 
gon. We estimate its D-perimeter from above using the estimates on A from the 
previous proof. 



ENI< E INI < (^ + o(\- 1 \o g \)) a 3 / \\ z \\dz<z. 

1 zt\D V7r ' Jd 



6 

INI < ( 

We need to estimate llvoll as well 



n 3 ' 2 , 



IKII < I -«ol = N < II -v \\ = HE^II ^ E INI ^ ™ 3/2 - 

1 

This shows that, indeed, L n <C n 3 ^ 2 . □ 

We mention that for a symmetric norm and for even n, the n shortest vectors 
can be chosen in pairs z, —z which is clearly optimal for L n . The case of odd n 
only causes only a minor difficulty 

Corollary 4.3 lim inf n~ 3 / 2 L n exists and equals a = a(D) > 0, say. 

5 Connection between L n and VP(tq) 

Lemma 5.1 Assume S G S with AreaS* = \, g{S) = 0. Then r G J- ' , the radial 
junction of S, is a feasible solution to VP(r ). Moreover, there is Q n G V n (for 
every n > 3) with 

7i r 7r r 3 (t} 

limn" 3/2 Per(X = -= / \\z\\dz = — = / —-4-dt. 

V6Js u " Zy/Gh r Q (t) 
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We remark that the last identity follows from a simple integral transformation. 



Proof. Feasibility of r is evident. We want to prove that for all e > (that we 
will suppose small enough without restricting the generality), there is Q n G V n , 
for every n > 3, with 



7T 

76 Js 



/ ||z||Gte — s < liminf n 3//2 Per Q n < limsupn 3 ^ 2 Per Q n < —= / ||z||dz + e. 
Js " V6 Js 



Since we would like to deal with sets whose perimeter can be defined and con- 
trolled, we introduce, for all m > 3, the m-gon approximation of S, whose ver- 
tices are r(— )u(— ) for k = 0,...,m — 1, recall that u{t) = (cost,sint). The 
sequence S m converges uniformly to S as m goes to infinity. Moreover, since 
g(S) = G int S, there are constants ci, c% > such that, for m large enough, 

ci-B c 5 m C c 2 B. 

We fix now m large enough so that the above condition is satisfied, as well as 

zdz < ce 



6 



TT 2 JS„ 



7T 

76 



7T 



Izlldz 7= 



\z\\dz 



< ce 



(Area S m — 1| < 



where c is a positive constant depending only on S that will be adjusted later. 

Now, there is a minimal A > (depending on m) so that |P fl \S m \ > n. Let 
L m denote the Euclidean perimeter of S. There are at most XL m lattice points 
on the boundary of \S m . Then, formula (13.61) applies and shows that 

|PnAS m | = ^AreaS' m + 0(A- 1 logA)) A 2 , 



implying A = n\Jn/ (6AreaS' m )(l + o(l)). 

We apply formula (13. 8p to \S with f(z) 
and f(z) = y where z = (x, y) to get 



or more precisely with f(z) 



x 



zePn\Sm 



6 



71- 



zdz + 0{\ 2 log A) 



Let P fl \S m = {zi, . . . , zi} (of course I > n) and define z = — J2[ %i. The 
previous equality implies that for n large enough ||^ || < 2ce\ 3 . The increasing 
slope construction applies to {^o, z%, . . . , z{\ and gives a convex lattice / or I + 1- 
gon T n . Note that T n has a special edge, the one parallel to, and having the 
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same direction as, Zq. All other edges of T n are short, shorter than C2G?2A <C n 1 ^ 2 
in D-norm. We claim now that, for a suitable choice of c (depending only on S), 
and e small enough, 

—j= I \\z\\dz — e < liminf n _3 ^ 2 PerT ra < limsupn~ 3//2 PerT n < —= / llzlldz + s 
V 6 J 3 v6Js 

We use again (13.81) this time with f(z) = \ \z\\ to get 

The claim follows since PerT n differs from Y^T I Nil D Y ||^o|| < 2ce\ 3 . 

Finally, let Q n be the convex hull of n consecutive vertices of T n , including 
the two endpoints of the special edge. Then Q n G V n and PerQ n < PerT™ 
and also, PerQ n is at least PerT™ minus the sum of the D-length of the missing 
edges, which is <C n as one can easily check. Thus |PerQ n — PerT n | <C n. The 
requirements on the constant c are now clear. □ 

We mention here that Lemma 15.11 implies Claim 14.21 by simply choosing any 
S G S with g(S) = and Area 5* = 1, for instance the Euclidean disk centred at 
the origin and having area 1. 



6 Bounding C n 

Our next target is to give bounds on the width and diameter of C n = convE^. 
Claim 6.1 The width of E n , w(E n ), satisfies w(E n ) 3> n 1//2 . 
Proof. Set w = w(E n ). Clearly, 

L n= INI > J2 M ^ M »W' 

where M n (w) is the sum of the lengths of the n shortest (in Euclidean norm) 
distinct vectors in Z 2 lying in a strip of width w. 

A simple yet technical computation, delayed to Appendix 1, shows that 
w < 772 1 / 2 (where 7 G (0,1/2]) implies M n {w) r?l 2 j^. This finishes the 
proof of Claim EIH because then n 3//2 ^> L n ^> M n {w) ^> n 3 / 2 /7 would lead to 
contradiction if 7 were too small. 

□ 



Claim 6.2 Assume the smallest Euclidean ball centred at and containing E, 
is RB. Then R < n 1 ' 2 . 
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Proof. Assume a is the farthest point (in Euclidean distance) from the origin 
in E n . Then |a| = R. Claim |4~21 implies that \a\ < L n <C r?l 2 . Since w(E n ) ^> 
n 1 / 2 by the previous claim, there is a point b £ E n whose distance from the line 
{x = ta : t G R} is > \w(E n ) > n 1 / 2 . 

The perimeter of the triangle A = conv{0,a,6} is \a\ + |6| + \a — b\ < 4|a| 
because |6| < \a\ and \a — b\ < \a\ + \b\ < 2\a\. Here Area A = \\a\h where h is 
the corresponding height of A. Since w(E n ) > n 1 / 2 , h ^> n 1 ^ 2 . 

Then by (13. 2p for large enough n, 



1 a , 6 A 1 A 
P n - A -Area -A 

' 2 ' 7T 2 2 



<3-2|a|log2|a| < /i|a|^P <AreaA-° g// 



n \/n 



implying that |P fl |A| > |Area A, again when n is large enough. 

Assume now that Area A > 16n. Then |P fl |A| > In. Since \E n \ < n, |A 
contains two distinct elements x, y G P\E n and, evidently, x + y G A. Then 
x,y,x + y G conv {0, a, 6, a + 6} contradicting Lemma 13.31 

Thus Area A = \\a\h < 16n, and so = |a| n 1 / 2 . □ 

We need one more fact about C n : 

Claim 6.3 Assume rB is the largest ball centred at and contained in C n . Then 
r 3> n 1 / 2 . 

Proof. Let a be the nearest point to on the boundary of C n . Thus r = \a\. 
Define E + = E n n {x G R 2 : ax > 0} and E~ = E n f] {x E R 2 : ax < 0}, and 
set f(x) = ax/\a\ which is just the component of x G R 2 in direction a. To have 
simpler notation we write f(X) = J2xex f( x ) when X G R 2 is a finite set. Since 
Y, z&En z = 0, f(E + ) + f(E~) = (because f(z) = when az = 0). We will show, 
however, that \a\ < 'jn 1 ^ 2 , for a suitably small 7 > 0, implies that 

f(E + ) + f(E-) < 0. (6.1) 

Define E + = {x G R 2 : < f(x) < ^n 1 / 2 } n RB with R < n 1 / 2 from 
Claim[6J2J The density principle tells now that f(E + ) < f(Pf)F + ) < /(Z 2 nF+) 
and the last sum can be estimated as follows. Let Q(z) be the unit cube centred 
at z. Again, AreaQ(z) n F + > 1/4 for all z G Z 2 n F + . This implies that, for 
large enough n, 

m := |Z 2 n F + \ < AreaF + /4 < i?|a| < 771. 

We use now (13.31) : 

/(z 2 nf + )-/ /W^ 



< -Rial. 



It is easy to see that f F+ f(z)dz < |a| 2 i? implying that /(Z 2 fl F + ) < |a| 2 i? -C 



7 2 n 3/ ' 2 . 
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Define F~ = {x E R 2 : > f(x) > -\^n 1/2 } fl RB where A > is chosen so 
that F~ contains exactly n — m — k lattice points. Here k is the number of lattice 
points on the line ax = so k < 2R + 1 n 1 ! 2 . Note that X^n 1 / 2 <C R since 
E n C RB consists of exactly n vectors. Choosing 7 small enough guarantees that 
m < O.ln which, in turn, guarantees that A > 1 and further, that |F~ fl Z 2 | > 
0.8n. The Euclidean perimeter of F~ is at most AR + X'yn 1 / 2 <C R and (13.11) 
shows that ||F~ fl Z 2 | — AreaF~| <C R. Clearly AreaF" <C RX'jn 1 ^ 2 , implying 
that 

0.8n < \F~ n Z 2 | < |l +0 ^ A7 ^i /2 ^J AreaF" < RX^n 1 ' 2 < A7n, 

which implies A7 ^> 1. 

The density principle says now that f(E~) < f(F~) (note that / is negative 
on F~ and E~), and f(F°) can be estimated using (13. 3ft : 



/(**-) - / /(*)cte 



< iT < n, 



because max{|/(x)| : x G F } < -R. Now /(z) is negative on F . It is easy to 
check that X 2/ y 2 nR <C — f F - f(z)dz <C X 2/ y 2 nR. So we have 



-f(F~)> [ -f{z)dz + 0{n)^> { -f{z)dz^>X 2 1 2 nR^> 
Jf- Jf~ 



n 3 ' 2 



This shows that ( 16. ID indeed holds if 7 > is chosen small enough because 
< /(Z 2 nF+) < 7 2 n 3 / 2 and -/(Z 2 HF-) >n 3 / 2 . □ 



Corollary 6.4 There are positive numbers r and R (depending only on D) such 
that for all n > 3 

rB C (AreaC n )~ 1/2 C n C RB. 



7 Almost all primitive points of C n are in E n 

We begin by stating a geometric lemma which is about a special kind of approx- 
imation. The technical proof is postponed to Appendix 2. 

Lemma 7.1 Assume K 6 K, is a convex polygon with rB C K C -RB. T/ien 
/or every 5 e (0, 0.02(r/i?) 2 ] t/iere are vertices V\,...,v m of K such that with 
Q = conv {r>i, . . . , v m } the following holds: 

• QCK C (l + AR 2 r- 2 6)Q, 

• for all i, the angle ZfjOfj+i is at least S. 
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Lemma 7.2 For every e > there is tiq = no(e,D) such that for all n > tiq, 



Proof. Let r n , resp. R n be the maximal, minimal radius such that r n B C 
C n C R n B. It follows from Claims [67TI and [6731 that R n jr n < c with a suitable 
positive constant depending only on D. Thus Lemma 17.11 can be applied with 
K = C n and 5 = e/(8c 2 ) (if e < 0.02/8 which we can clearly assume). We get a 
polygon Q = conv {vi, . . . , v m } satisfying C n C (1 + e/2)Q. 

Assume, contrary to the statement of the lemma, that there is an x G (1 — 
e)C n nP\£ n . One of the cones pos {vi, fj+i} contains x, say in the cone W : = 
pos{vi,v 2 }. Define A = conv {0, Vi, v 2 }. Thus A C C n n W C (1 + e/2)A. 
As x e (1 - £)C n n W, ui + w 2 - ac G W \ (1 + e)A. The triangle A* = 
((t>i + v 2 — x) — W) \ (1 + s/2)A is disjoint from C n . We claim that it contains 
a primitive point y. This will finish the proof since then x,y,x + y all lie in the 
parallelogram with vertices 0, V\, v 2 , V\ + v 2 contradicting Lemma [3.31 

We prove the claim by using (13. 2ft : Area A* ^> e 3 n because its angle at 
v i + v 2 — x is at least S, and the neighbouring sides are of length at least £|vi|/2 
and e\v 2 \/2 and \vi\, \v 2 \ ^> n 1 ^ 2 . Further, its perimeter is at most |t>i| + \v 2 \ + 

— v 2 \ n 1 ! 2 . Thus 



(1 - e)C n nPc£, 



A* n P 



6 



-Area A* ^ (log n)n 



1/2 



7T 



Here Area A* is of order e 3 n and the error term is of order (log n)n 
e fixed, A* contains a primitive vector if n is large enough. 



I 2 . Since 



□ 








(l-£>i «i (l+e/2)«i 



Figure 2. The proof of Lemma 7.2 



8 Proof of Theorem E3 



In this section we prove Theorem 12.21 apart from the uniqueness of C and r which 
will be shown in the next section. 



I. Bar any, N. Enriquez - Jarmk's convex lattice n-gon 



13 



The Blaschke selection theorem and Corollary 16.41 imply that every subse- 
quence of (AreaC n )~ 1 / 2 C n contains a convergent (in Hausdorff metric) subse- 
quence. Corollary 14.31 guarantees then the existence of positive integers ri\ < 
n 2 < . . . such that limn fc 3//2 L„ fe = a and limdist ((AreaC nk )~ 1 ^ 2 C rik ,C) = for 



some convex body C G JC. Define X k = ^AreaC nk and set, for simpler writing, 

C k = X k l C nk . It is evident that rB C C C RB, showing that, for every 5 > 0, 
{l-5)C C C k C (1 + 5)C for all large enough k. Since n k = 4,AreaC„ fc (l + o(l)), 
A fc = ^v^(l + o(l)). 

It follows immediately that AreaC = 1. We show next that J c zdz = 0. 
For this it suffices to prove that J c f(z)dz = in the case when / is the linear 
function f(z) = x and f(z) = y where z = (x,y). Choose e > and then, using 
Lemma [T7T| k Q so large that, for k > k Q , 

(1 - e/2)C nk n P C E nk c C nk n P. 

It follows now that there is a k\ so that for all k > k\ 

(1 - e)\ k C n P C E nk c (1 + e)\ k C n P. (8.2) 

Using the notation /(X) = X^ 6 x /l- 2 ) when X C R 2 is finite, we have f{E nk ) = 
0. Next, 

|/(PnA fe C)| = |/(PnA fc C)-/(£ n J| 

< |/(Pn[(l + e)A fc C\(l-e)A fe C])| 
<C eA^ max{/(z) : z G A^C} <C 

On the other hand, by (13.81) . 

|/(P n X k C)\ = ^X\ J c f(z)dz (l + G>(A^ log A*)) 

as one can check easily. So if J c f(z)dz ^ 0, then /(P n A^C) is of order n k . 
But as we have just shown, |/(P fl X k C)\ <C en k . So indeed, J c f(z)dz = 0, or, 
in other words, g{C) = 0. 

An almost identical proof, this time with the 1-homogeneous function f{z) = 
\\z\\ gives 

I \\x\\dx = a(D). 
We only give a sketch: Equation ( 18.21) shows that 



E INI- E IN 



<C en k . 
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Here E 2e pnE„ fc INI = L n k and so \imn k 3/2 E ze pnA fc C INI = 01(D). The estimate 
(13. 4p says now that 



H^ll ( 

jePnAfcC 71 JXkC 



< n fc logn fc , 



and J c ||x||<ix = a(D) follows. 

Lemma 15.11 applies now because g{C) = and AreaC = 1. So there is a 
sequence Q n £ V n with limn _3 ' /2 Per Q ra = a(D). Then L ra < PerQ n implies that 
limn-3/2^ _ a (D). □ 



9 The variational problem 

Next we turn to uniqueness. As first step we treat a special case. 

Lemma 9.1 Let r$ be the radial function of D e /C with g(D) = 0. T/ien ro 
t/te unique solution to VP{tq). 

Proof. We consider the variational problem which ignores the constraints 
about the center of gravity : 



r Z7T 

minimize / r 3 (t)/r (t)dt 
Jo 

subject to / r 2 (t)dt = 2 
Jo 



From Holder's inequality : 

2-7T / p2tt r 3 

\ 2 / 3 / r2n \ 1/3 

c r2< ~ U ^) (/„ r 

which is an equality if and only if r and r are proportional. In our case J 21T r 2 = 
J 2f rg = 2 and so r = ro . □ 

We now use the previous lemma to treat the general case: 

Lemma 9.2 There exists a unique solution r G T to 'problem VP{tq). This 
solution is equal to 

( a . N_1 

r = ho cos t + c sin t 

Vr 

where a > 0, 6, c are i/ie unique real numbers which make the function r satisfy 
the three constraints of VP(r ) . 
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Proof. We prove in Appendix 3 that every optimal solution r G T c to VP(ro) 
is of the form r(t) = (^ L + 6cost + csint) _1 with suitable constants a, b, c G P. We 
have shown that the radial function, r(t), of C from Theorem 12.21 is an optimal 
solution to VP(r ). As C is convex, r(t) is equal to + boost + csint) -1 . 
According to the previous Lemma, the unique solution to the variational problem 
VP(r) is r. 

Consider now another optimal solution, r*, to VP(r ). It is clear that r* is a 
feasible solution to VP(r) and that 



Further, 



o r Jo r 



a / — = / r* 3 h 6 cos t + c sin t) = — , 

Jo r Jo \r J Jo r 



and, in the same way, 

r 2n r 3 f 2n , / a \ / >27r r 3 

, / — = / r h boost + csint = / — . 

Jo r 7o Vr J Jo r 



So r*, too, is an optimal solution to VP(r). By the Lemma, r = r*, and a > 
follows as well. □ 

Remark: After reading this proof, one easily understands that r(t) is the 
unique solution to the variational problem in a class of functions larger than T . 



10 Proof of Theorem 

This is fairly simple once we know that C is unique. Let u(t) = (cost, sin t) be 
the unit vector in direction t G [0, 27r]. When a minimizer P n is translated as 
Theorem 12.31 specifies, the sum of the edges of P n having direction between u(0) 
and u(t) is very close to the sum of the primitive vectors having direction between 
u(0) and u(t) in C n . The latter, divided by n 3//2 is very close to P(t) = J c ^ zdz 
where C(t) is the set of vectors in C with direction between u(0) and u(t). The 
curve P{t) is closed (because g(C) = 0) and convex (this has been shown in [2]), 
so it is the boundary of a convex set P. The simple and straightforward checking 
of 

limdist (n~ 3/2 P n ,P) = 

is left to the reader. We remark that the convexity of P(t) follows also from the 
fact that the boundary of P n , after suitable rescaling, tends to P(t). □ 

The same construction C — ► P with P(t) = j C i t \ zdz is used, with a similar 
purpose, in [2]. Further properties of the construction are also established there. 
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11 An example 



We concentrate now on the cases when the solution is constant which correspond 
to the case when the limit shape of the polygon is a circle. 

Lemma 11.1 The solution is constant if and only if 1/tq is of the form a + 
bcosO + csin6> ; or, in other words, when r is the radial function of an ellipse 
having its focus point at the origin. 

Proof. Suppose the solution is constant, the form of ro is then directly 
derived from Lemma 9.2. Conversely, if ^ is of the form a + bcosO + csin#, the 
solution is then also of the form - = a' + b' cos 6 + d sin 6. This says that it is the 
radial function of an ellipse having its focus point at the origin. We conclude by 
observing that the only ellipses whose centre of gravity is at the same time their 
focus point, are circles. □ 

12 Appendix 1 

Lemma 12.1 Let M n (w) be the sum of the lengths of the n shortest (in Euclidean 
norm) distinct vectors in Z 2 lying in a strip of width w, centred at the origin. 
Suppose 7 G (0, 1/2], then w < 771 1 / 2 implies M n (w) ^> n 3 ^ 2 /^. 

Proof. It is clear that this set of vectors is just the set of lattice points 
contained in A := dB n T where T is a strip of width w, centred at the origin, 
and d is a suitable radius making A fl Z 2 have exactly n elements (ties broken 
arbitrarily). Let (p denote the angle that the strip T makes with the x-axis of 
R 2 . We may assume by symmetry that (p G [0, 7r/4]. 

Observe first that d > ^/n/2 since otherwise the disk dB would contain fewer 
than n lattice points. Let Q(z) denote the unit square centred at z G R 2 and let 
£k be the line with equation x = k (k is an integer). Clearly, £k intersects S in 
a segment of length wcosip, and so it fl Z 2 contains at least [w/ cosy?J and at 
most [w/ cosy?J + 1 lattice points from S. 

Assume first that wj cosy? > 1. As is easy to see, Area A fl Q(z) is at least 
1/4 for z G A n Z 2 . Hence, Area A > n/4. Since Area A < 2dw, d > n/(4w) 
follows. 

For simpler notation write u = (dcosip)/2. For the lines £k with k G [u, 2u — 
w/2], £ k fl A contains at least [w/ cos(/?J lattice points. Since w < u, there are 
at least \2u — w/2\ — \u\ ^> u such lines. All of them have distance at least 
(d — w)/2 ^> d from the origin. Consequently, using the bounds w < ^n}l 2 and 
d > n 1 / 2 /2 generously, 




I. Bar any, N. Enriquez - Jarmk's convex lattice n-gon 



17 



Assume next that wj cosif < 1. There are at most six z G A D Z 2 such that 
Q(z) intersects the boundary of dB. For the other z G AnZ 2 , Q(z) intersects the 
boundary of A in one or two line segments, whose total length is between 1/ cosip 
and 2/ cos^j. For distinct lattice points in A fl Z 2 the corresponding segments do 
not overlap. This implies that 

^<4rf< 2(n ~ 6) . 

COS (fi COS </J 

Each line with < 2u/3 contains at most one lattice point from A. The 
remaining points from A fl Z 2 , and there are at least n — 2|_2w/3j — 1 of them, 
are at distance ~ — 1 from the origin. Hence, we see 

M„ W > g - l) (n - L 2^j - l) » r?. 

□ 

13 Appendix 2 

We start the proof of Lemma 17.11 with the following Claim. 

Claim 13.1 Suppose a,b,c,d are vertices of K (in anticlockwise order), [a,b] 
and [c, d] are edges of K , and IbOc < 35. Let x be the intersection point of the 
lines through a, b and c, d, and let y be the intersection point of the lines through 
0, x and a, c. Then \x — y\ < A5(R/r) 2 \y\. 

Proof. The condition rB C K C RB implies that /3 = lOxb = /.Oba — Lxba > 
arcsin r/R — 36 since 

sm lOba = 

\b\ 

{£a,b being the line through a and b) \b\ < R , d(0, l a ^) > r by assumption, so 
that sin lOba > r/R, see Figure 3. 

Further Lxyc = lOxa — 1x0b > (3. The sine theorem in the triangle x, y, c 
shows that 

\x — y\ sin Lcxy 
\x — c\ smLcyx'' 
and similarly, the sine theorem in the triangle x, 0, c shows that 

| a; — c\ sin IcQx 
\x\ sin lOcx 

Multiplying them gives 

| a; — y\ sin Lcxy sin IcOx sin 35 
\x\ sin Lcyx sin LOcx (r/R) sin (3 
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Next, since \y\ 



\x\ 



\x 



y\, we have 





X 




\x 




\x - y\ 



\y-y\ 

\x\ 



< 



sin 3<5 
(r/R) sin/3 



We use this inequality next in the form 
\x - y\ 



< 



\y\ 



sin 35 


X 




(r/R) sin/3 


X 




\x - y\ 



< S ^ n ^ < 45 ( ^ 

(r/R) sin (3 — sin 35 \r 



where we only have to check the validity of the last inequality. This is a matter of 
direct computation using that sin/3 > sin(arcsin(r/i?)— 35) > (r/R) cos 35— sin 35 
and the assumption that 5 < 0.02(r/R) 2 implying, in particular, that 5 < 0.02. 
What is to be checked now is that 



tan 35 



1 + 45 



r 
R 



< 45. 



Here 5(R/rf < 0.02 and so the expression in the square bracket is at most 1.16 
and the inequality follows. We omit the details. □ 




o 



Figure 3. The proof of Claim 13.1 

The Proof of Lemma 17.11 is an algorithm that constructs the vertex set V of 
Q. We start with V = 0. We call the edge [a, b] of K special if laOb > 5. Let 
W be a cone with apex at and angle 5. It follows that if W is disjoint from all 
special edges, then it contains a vertex of K. 

Case 1. Let [aj, bi], [a 2 , b 2 ], ■ ■ ■ , [afc, be consecutive special edges in an- 
ticlockwise order so that Z6j0a i+1 < 35 for alH = 1, . . . , k — 1 (or up to k if 
ZfefcOa! < 35). We call this a maximal chain of consecutive special edges if there 
is no special edge [a, b] with lWa\ < 35 or IbkOa < 35. 
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For such a maximal chain we put the vertices a±, . . . ,af.,bf. (or a%, . . . , if 
Z&fcOai < 35) into V, and we do so for all such maximal chains. 

Case 2. Let [ai,6i] and [02,62] be consecutive special edges with vertices 
01,61,02,62 in anticlockwise order so that 7 := Z6i0a2 > 35. Then we choose 
5' G [5,35] so that 7/0 7 is an odd integer, say 2h + 1. This is always possible 
since there is an odd integer between 7/(35) and 7/5 because their difference is 



Subdivide now the cone pos {61, 02} into 2h + 1 subcones, each of angle 5' and 
pick a vertex u\, . . . ,Uh from every second subcone. Finally, put 61, U\, . . . , Uh, a 2 
into V. 

If there are only two special edges [ai, 61] and [02, 62], then one has to do the 
same construction between edges [02,62] and [01,61] as well. If there is only one 
special edge, then the construction is carried out from 61 to ai as if one had two 
special edges [01,61] and [61, 01]. 

Finally, if there are no special edges, then we chose a 5' G [5, 25} so that 
2tt/5' is an even integer, 2h, say. This is evidently possible. Subdivide the plane 
into cones of angle 5' (with apex at 0) and choose a vertex u\, . . . ,Uh from every 
second cone, and set V = {ui, . . . , Uh}. 

The algorithm is finished. By construction ZfjOfj+i > 5: for the angle at 
0. Finally we check condition K C (1 + 45(R/r) 2 )Q. Let t>j, v i+ i, v i+2 , ^1+3 be 
four consecutive vertices of Q in anticlockwise order. Rename these points as 
o, 6, c, d as in the Claim. Then K R pos (6, c) \ Q is contained in the triangle 
6, c, x from the Claim. Now y G Q because y lies on the segment [o,c], and so 
x G (1 + 45(R/r) 2 )Q according to the Claim. So the triangle 6, c, x is contained 



14 Appendix 3 

It happens that standard theorems of the Calculus of Variations (see for instance 
[6]) are stated in a C 1 setting, and suppose also that the function r$ involved in 
the problem is C 1 . Since these conditions are not satisfied in our problem, we 
have to elaborate the following statement: 

Lemma 14.1 All the solutions r G T c satisfying problem VP(r ) are of the form 



where a, b, c are real numbers which make the function r satisfy the three con- 
straints of problem M.ty) . 

Proof. Consider r an optimal solution in JF C . Let h be a function on [0, 2ir] 
such that the perturbed function r e := r + eh remains in T c for e in a neighbour- 
hood of (notice that all the twice different iable functions are convenient). This 



1 /5- 1 /(35)= 2 1 /(35)>2. 



(l + 45(R/rf)Q. 



□ 



r = h 6 cos t + c sin t 

Vr 
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perturbation won't be feasible in general. We want to modify it in order to make 
it feasible. That is what we do in the two first steps. 

Step 1. We translate the set defined by the function r £ in order to get a 
centred set defined by a new radial function f £ we evaluate up to some o(e). 

In the following, the notation o(e) stands for some family of functions, which 
may be constant, indexed by e, such that both ^ converges to as e goes to 0, 
and ^ is dominated. 

e 

The coordinates of the centre of gravity of the set defined by r £ are 



( f 2lT s f 2w s 

\J r'l cos tdt, J r £ smtdt 



/•27T c-lft \ 

J 3r 2 h cos tdt, J 3r 2 hsmtdtj + o(e) 



Recall that u(t) = (cost, sin t). Define the numbers r h , 9h by setting r h u(9 h ) '■ = 
(J 2n 3r 2 h cos tdt, J 2w 3r 2 h sin tdt) . 

For a given 6, the polar coordinates of r £ u(9) — er h u(9 h ) + o(e) are given by 

r(9)=rM-er h cos(e h -e) + o(e) 
Hence, f can be expressed as a function of 9 as follows: 



r(9) = r £ 



9 + £r h Sm ^ j ^ + o(e)^ " er h cos (9 h -9 + o(e)) + o(e) 



Using now the almost everywhere differentiability of r (and therefore of r £ ) which 
is inherited from convexity, we obtain that, almost everywhere, 



f{9) = r(9)+e 



h{9) + r h sm(9 h - 9)-{9) - r h cos{9 h - 



+ o(e) 



Note that the domination of in the last step is due to the fact that the left 
and right derivatives of r are bounded on [0, 2tt]. 

Step 2. We obtain a completely feasible function r* , by normalizing f by the 
area of the set defined by f, which is the same as the area of the set defined by 
r £ , since the two sets are obtained one from the other by a translation. 

Define, 



r f (9) 



f{9) 



f(9) 



The function r?(9) can be written as r{9) times the function 



(9) ^-l(J Q 2n rh) + o(e) 



l + e 



r(9) 



+ r h sin(9 h - 



r'(9) 



- r h 



cos(9 h -9) 1 r*K 



r{9) 



1 r Z7T 
-2</„ rfc > 



+ o(e) 
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Step 3. Now, we test the optimality of the function r by considering the 
functional applied to the feasible perturbation and writing the integral J 2n 

r 27r rf 



as P lus 



3e 



Jo r n 



r(9) 



r 2 (9) 



r(9) 



271 



rh) 



+ o(e) 



When developing the sine and cosine in the above bracket and performing the 
integration on 9 (and keeping in mind that r^ and 9h are constants that don't 
depend on 9 !) we deduce that, if r is optimal, there exist real constants A, B 
and C such that for all twice differentiable function h, 



p2tt rf^Yi f2n 

/ + Ar h sin 6 h + Br h cos 9 h + C / rh = 

Jo r Jo 



Recall that (r h cos9 h , r h sm9 h ) = (J 2w 3r 2 h cos tdt, J 2w 3r 2 h sin tdt). 
Therefore, for all twice differentiable functions h 



2tt / I 

hr 2 h 3 A cos i 



3fisin^ + - 

Jo \r r 

which implies that the bracket inside the integral is 0. 







□ 
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